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Abstract. The present-day state of the Milky Way disk can tell us much about the 
history of our Galaxy and provide insights into its formation. We have constructed a 
high-precision catalogue of disk stars using data from the Sloan Digital Sky Survey 
(SDSS) and use these stars to probe the heating history as well as investigating the 
detailed phase-space distribution. We also show how this sample can be used to probe 
the global properties of the Milky Way disk, employing the Jeans equations to provide 
a simple model of the potential close to the disk. Our model is in excellent agreement 
with others in the literature and provides an indication that the disk, rather than the halo, 
dominates the circular speed at the solar neighborhood. The work presented in these 
proceedings has b een published as "Slicing and dicing the Milky Way disc in SDSS" 
(ISmith et al]l20l2h . 



1. Introduction and sample construction 

Large surveys such as SDSS are beginning to enable detailed studies to be made of the 
properties of the Milky Way. It is now possible to construct high-precision catalogues of 
stars with full 3D kinematics and from these analyse the structure and evolution of our 
Galaxy. In the following p roceedings I sum marise our study of the disk of the Milky 
Way (presented in full in ISmith et al.ll2012h . highlighting two of the most important 



results. 

We use data from the 7th public data release of the SDSS, concentrating on the 250 
square degree Stripe 82 region for which it is possible to obtain high precision photome- 
try and proper motions (Bramich et alj|2008). This catalo gue has been exploi ted in a va- 



riety o f works, covering white dwarfs (Vidrih et alj|2007l). hal o kinematics (Smith ^t al.l 
2009al lbh and RR Lyrae in the outer halo <|Watkins et al. 200^). 



We take all dwarf stars with good quality spectra dLee et al. 20081) . estimating dis- 



tances using a slightly modified version of the llvezic et al.1 (120081) photometr ic parallax 



relati on. Our modification concerns the turn-off correction (equation A6 of Ivezic et al 



2008), which we argue is not applicable to disk stars (see the Appendix of ISmith et al 



20121 ). We then calculate the three-dimensional velocities for our sample, propagating 
all uncertainties, and obtain final errors of around 20 to 30 kms" 1 for each component 
of the velocity. 

As we are interested in the kinematics of the disk as a function of height from the 
plane, we restrict ourselves to 7 < R < 9. We also split the data into three ranges in 
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metallicity (-1.5 < [Fe/H] < -0.8,-0.8 < [Fe/H] < -0.5 and -0.5 < [Fe/H]), and 
then for each metallicity bin we further divide the data into four ranges in z out to a 
maximum distance of 2 kpc. The stars are equally divided between the four distance 
bins, resulting in in around 500 to 800 stars per bin from our full sample of 7280 stars. 
This binning is necessary in order to avoid having to model the SDSS spectroscopic 
selection function. 

Once we have divided up our stars, we then proceed to fit the distribution of kine- 
matics using maximum likelihood methods (incorporating the uncertainties on the in- 
dividual velocities). Both vr and v z are fit using Gaussian distributions, but since it is 
well known that the distribution of v^ for the disk is highl y skewed and non-Gaussia n 
we adopt an asymmetric functional form (equation (15) of lCudd eford & Bin ney|[T994 ). 
When carrying out this analysis it is crucial to remove contamination from halo stars, 
particularly for the most metal-poor bins. We do this by including a non-rotating Gaus- 
sian component into the velocity fits, where the level of contamination (and hence am- 
plitude of the Gaussian) is assumed to be twice the number o f counter-rotating st ars. 
The halo dispersions are fixed according to the values found in lSmith et al.1 (12009a ). 



2. The rotation lag of the disk and the asymmetric drift 



A comprehensive analysis of the results are presented in lSmith et all (120121) . Here we 
concentrate on the most important properties, starting with rotation lag. We are able 
to trace the lag (Fig. [D to a couple of kpc outside the plane and find that it follows 
the well-known asymmetric drift relation. The hotter (metal-poor) stars exhibit greater 
lag than the colder (metal-rich) populations, resulting in a clear correlation between 
lag and metallicity. The gradient of the lag with respect to |z| varies from around 15 
to 40 kms _1 kpc _1 , depending on metallicity. Fig. [T]also shows how the lag is cor- 
related with cr~ . For the solar-neighbourhood it is know that this is a linear relation 
dDehnen & Binne3[l998ab, but as we can see from our figure this is no longer true once 



we move beyond 0.5 kpc. However, it can be seen that there is still a relatively tight 
correlation with <r|, which is independent of metallicity. 



3. Constraining the gravitational potential of the disk 



As has been shown by many authors, most notably in the seminal work of lKuijken & Gilmore 



( 1991), the vertical velocity dispersion profile can be used to constrain the gravitational 
potential of the disk. We show the results derived from our data in Fig. |2l where we 
have included an additional sola r-neighbourhood metal -rich data point derived from 
the Geneva-Copenhagen survey dNordstrom et all [20041) . We take a simplified two- 
parameter model for the potential, consisting of an infinite razor-thin sheet embedded 
in a constant background. Our best-fit model estimates that the mass density of these 
two components are 32.5 M Q pc~ 2 and 0.015 M Q pc~ 3 , respectively. 

As can be seen from Fig. |2l our measurement of the potential is in good agree- 
ment with existing models, especial ly Model 1 of Dehn en & Binnevl ( 1998bl) . This 



is the least halo-dominated model of iDehnen & Binnevl (Il998bl) . with the disk (which 



has a relatively short scale-length) dominating the cir cular speed at the solar neig hbor 



hood. Interested readers should consult section 2.7 of Binney & Tremain e (2008) for a 
detailed comparison of these models. 
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Figure 1. The rotational lag, plotted against z (left) and against the radial ve- 
locity dispersion (right). The triangles, squares and crosses correspond to metal- 
poor, intermediate-metallicity and metal-rich populations, respectively. The dot- 
ted line in the right pane l corresponds to the solar-neighborhood relation from 
iDehnen & BinnevT(ll998al) and the solid line denotes an empirical fit wi th the lag 
equal to 0.0149o-2 + 1.21 x 10" 6 cr£ • Figure taken from lSmith etal] (120121) . 



Note that this background mass density is close to t he 0.014 M n pc~ 3 predicted us- 
ing isothermal spherical halo models (equation 4.279 of iBinney & Tremainell2008h . If 
we assume our background mass represents the dark halo, it corresponds to a local dark 
matter density of 0.57 GeV cm -3 , which is noticeably larger t han the canonical value 
of 0.30 GeV cm~ 3 typically assumed (e.g. lJungman et aJll996l) . Our analysis adds still 
more weight to the argument the local halo density may be substantially underestimated 
by this generally accepted value. Perhaps more robust than the local mass density is the 
surface mass density. By integrating our mass distribution we obtain a total surface 
mass density of Si.i kpc - 66 M pc~ 2 , w hich agrees well with the classical value of 
71 ± 6 M pc~ 2 from Kuijken & Gilmord (Il99ll) . If we integrate beyond 1.1 kpc, we 
find I 2 k P c - 94 M Q pc -2 and Z4k pc = 155 M Q pc~ 2 . 



4. Conclusion 



Our results (presented in full in iSmith et all 120121) can be used to address the global 



properties of the disk and also the heating processes that have shaped its evolution. In 
order to fully understand the nature of disk heating and to disentangle the contribu- 
tions from various mechanisms, one needs to go beyond the work presented here. The 
most crucial improvement will be the ability to make direct estimates for stellar ages, 
rather than relying on correlations with metallicity. One aspect that will help us in this 
effort is by folding in measurements of alpha-elem ent abundances that are now being 
determined routinely for vast numbers of stars (e.g. lLee et al.ll201 lb . 

Question: (N icolas Martin) C an you develop your model to test for the presence 
of a dark disk (e.g. Re ad et alj|20 08)? Reply: Due to the limitations of the data that we 
have used there are a number of uncertainties in our modelling. As a consequence of 
these limitations, I think it would be difficult to include extra parameters. 
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Figure 2. The results of our modeling of the potential of the Galactic disk. The 
left panel shows the vertical velocity dispersion profiles (metal-poor - top; metal- 
rich - bottom), along with the corresponding profiles for the potential (solid lines). 
The right panel shows the potential resulting from this model. For the purposes of 
comparison we have included models for the potential taken fro mlDehnen & Binnevl 
(1998 - Model 1, dotted; Model 4, dashed). Figure adapted from lSmifh et al I (l2012h . 
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